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EVALUATION OF THE PRESENT THEORETICAL BASIS FOR 
DETERMINATION OF PLANETARY SURFACE PROPERTIES 


BY EARTH-BASED RADAR 

Leo D. Staton 
Langley Research Center 

SUMMARY 

Spaceflight programs such as the planned Viking landing on Mars often require the 
determination of planetary surface slopes and surface dielectric constants by Earth- 
based methods. Heavy reliance is often placed on radar backscattering data for estima- 
tion of these surface properties. The present report, which is primarily directed toward 
users of such estimates rather than toward radar workers, is an assessment of the basic 
theory by which the raw radar data are interpreted. Certain serious difficulties and inter- 
nal inconsistencies in the available theoretical formulas are reemphasized. These diffi- 
culties include questions of the fundamental applicability of Gaussian height statistics to 
planetary terrain, questions of the validity of certain assumptions necessary for applying 
results for perfectly conducting surfaces to dielectric surfaces, and questions surround- 
ing a widely used rough-surface autocorrelation function that leads to several absurd and 
self -contradictory conclusions. The discussion therefore brings into serious question the 
reliability of the presently available results for these surface properties as obtained by 
Earth -based radar methods, and it points out that the theoretical difficulties are not likely 
to be remedied in the near future. 


INTRODUCTION 

The interpretation of electromagnetic signals scattered from rough surfaces has 
important applications in both terrestrial and astronomical studies. This report concen- 
trates on certain aspects of the astronomical area, particularly on methods of inferring 
surface electrical and scattering properties from data obtained by Earth-based radar 
reflections from celestial bodies (for example, from Mars). Since they can reveal infor- 
mation about characteristic planetary surface slopes and surface particle composition 
and/or density, such inferences are important for the planning of planetary soft -landing 
missions such as the NASA Viking exploration of Mars in 1976. The principal aim of this 
paper, therefore, is to review and to criticize the basic theory by which radar astronomi- 
cal data have been interpreted in pursuit of this practical goal. Since rough-surface 



scattering theory is still a highly active field of research, it is not surprising that most 
existing inferences of planetary surface properties, such as detailed surface slopes and 
dielectric constants, rest on somewhat unsatisfactory grounds. Of course, this does not 
demean the spectacular successes of radar astronomy over the last 20 years in other 
areas of lunar and planetary studies, such as large-scale topographic mapping, planetary 
rotation-rate determination, and interplanetary distance measurements. 

In the present report the fundamental statistical notions and the electromagnetic 
scattering approximations which are the basis for most inferences of surface properties 
from radar astronomical data are discussed, and certain inconsistencies in this work are 
pointed out. Some of these inconsistencies have been discussed in the literature of 
rough -surface scattering but have not yet been taken into account in the planetary radar, 
work. 


SYMBOLS 

A covariance matrix defined by equation (37) 

C quantity defined as coefficient of sin^ <p in equation (47) 

D upper radial cutoff value associated with equation (44) 

E magnitude of electric field 

e magnitude of unit polarization vector 

F integrand function of equation (9); also, function defined by equation (38) 

f physical quantity associated with a statistically described surface 

H magnitude of magnetic field 

h local height of rough surface 

I unit dyad 

Jo zeroth-order Bessel function of the first kind 

j imaginary unit, \[^1 
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wave number, 2n/X 

total length and total width of scattering surface 
surface correlation length 
cofactor of covariance matrix 

total number of sample functions in a statistical ensemble; also, total number 
of specular points on a rough surface 

number of sample functions in ensemble having height between h and 
h + dh; also, magnitude of unit normal vector 

angular power scattering function at angle <p 

probability density 

reflection coefficient for plane -wave incidence on plane surface 
magitude of difference between incident and scattered propagation vectors 
distance from radar antenna to planetary center 
distance from point on scattering surface to field point 
area of scattering surface 

slope of surface in x -direction, ^ 

distance between point with coordinate x and a given stationary point on 
rough surface 

coordinates with respect to mean surface 
local slope angle 
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y 


angle between local surface normal and overall mean surface normal 


A operator indicating small increment 

e dielectric constant 

\ wavelength 

p magnetic permeability 

p radius of curvature; also, autocorrelation function of rough surface 

2 total radar cross section, length squared 

cr standard deviation; also, radar cross section per unit width, length; and radar 

cross section per unit area, dimensionless 

r time 

<p phase angle of individual scattered wave in equation (21) 

<P planetary great-circle angle from sub -Earth point, equal to angle of incidence 

with respect to mean surface 

Green's function of equation (6); also, local angle of incidence 
w angular frequency 

Superscripts: 

sc scattered field 

* complex conjugate 

' quantity evaluated on scattering surface; also, first derivative 

" second derivative 
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Subscripts: 


a,b two principal directions of curvature on a surface 

f property associated with the physical quantity f 

h horizontal polarization direction 

i general member of a collection; also, incident field 

j,k general member of a collection 

m,p specular point, summation index 

o total field incident on planet 

R received signal 

s scattered field 

sp quantity evaluated at specular point 

v vertical polarization direction 

x vector component in x-direction; also, derivative with respect to x in equa- 

tions (24) to (28) 

y vector component in y-direction; also, derivative with respect to y in equa- 

tions (24) to (28) 

z vector component in z -direction; also, general polarization direction 

1,2 dummy variables of integration 

Mathematical notation: 

- (over symbol) average value over a particular sample surface of an ensemble 
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(over symbol) vector quantity 


( ) expectation value computed over an ensemble of sample surfaces 

~ (over symbol) unit vector 

( ; ) simple probability density 

( | ) conditional propability density 

(||) j oint probability density 

V vector gradient operator 

1 1 absolute value of a quantity 

STATISTICAL NATURE OF THE PROBLEM 
OF ROUGH-SURFACE SCATTERING 

The exact calculation of the fields scattered by an object immersed in a specified 
incident electromagnetic field is an intractable problem for all but the simplest geomet- 
rical shapes, and solutions for even these simple shapes often require extensive numeri- 
cal computations, the complexity of which depends upon the ratio of the incident wave- 
length to the characteristic size of the scattering object. Thus, for scattering objects as 
complicated as rough planetary surface terrain, no exact solution is possible and the 
problem becomes a search for an appropriate approximation scheme. Furthermore, if 
the basic problem is to deduce the properties of the scattering object by measurements 
on the scattered fields, the difficulty is manifestly compounded in that no unique solution 
exists. That is, it is impossible, in principle, to define the scattering object solely from 
measurements of the scattered field taken at any finite number of points in space. For 
any progress at all, then, plausible models of the scattering object must be constructed, 
with the field measurements being used to fix the values of whatever parameters the 
model incorporates. For the study of scatterers having the qualities of rough planetary 
terrain, for which it is impossible to ascribe, a priori, any particular shape, it is clear 
that the model must be a statistical one, the results of which must be interpreted only in a 
statistical sense. 
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Statistical analyses of rough-surface scattering have been developed by many 
authors, most of them using the techniques of S. O. Rice (ref. 1). These techniques are 
very closely related to the methods used in statistical communication theory (ref. 2). It 
is only necessary to replace the time variable by a spatial variable and to perform an 
obvious extension from one to two independent variables in order to adapt the formalism 
of reference 2 to the description of the statistical quantity of interest for rough-surface 
scattering problems. This quantity is the variable terrain height defined over some 
datum surface such that this height is a random function of position on the surface. The 
following elementary ideas are reviewed in this section for reference. 

The term random function means a quantity, in this case the height, which takes on 
variable values for its various argument coordinates x and y, such that at a given point 
the quantity takes on a particular value within the interval (h, h+dh) according to a prob- 
ability p(h;x,y) dh. The situation is thus describable by a (sufficiently) large collection 
or ensemble, say N in number, of the physically possible realizations of specific func- 
tions Iq(x,y) (for i = 1, 2, . . ., N) such that at every point (x,y), each physically possi- 
ble value of h, to within some arbitrarily small tolerance dh, is taken on by some one 
or more, say n(h,x,y), of the lq(x,y). Then the probability density p(h;x,y) is related 
to this ensemble by the relation 

p(h ; x,y)dh=2<!!i>4 (1) 


The totality of information in the ensemble of specific functions tq(x,y) constitutes a 
random process. All this information may be represented by the infinite set of joint 
probability densities pj(hi,h 2 ». . .,hj||x 1 ,y 1 ;x 2 ,y2;* • defined by 

, „ \ n(h 1 ,x 1 ,y 1 ;h 2 ,x 2 ,y 2 ;. . .;h.,x 1} 

Pj ( h l, h 2,* • •>hj||x 1 ,y 1 ;x 2 ,y2;. . -;xj,yj) dh^ . . . dhj = jj J - J - 


y j) 


( 2 ) 


where j = 1, 2, 3, . . ., °°, and n(h 1 ,x 1 ,y 1 ;h 2 ,x 2 ,y2;. • .;hj,Xj,yj) is the number of 
functions in the ensemble which take on the specific set of values h^ at the points (x^y^) 
within the arbitrarily small tolerances dh^, and where k takes on the values 
1, 2, . . ., j. 

With these notions, physical experiments such as the following can be analyzed. 
Suppose that some physical quantity, say f(h,x,y), which is a function of height and posi- 
tion on a randomly rough surface, must be measured. Further, suppose that a collection 
of surfaces which have been produced by the same physical process is available for study, 
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each surface differing in the detailed height variations upon it. The fact that the surfaces 
were produced by the same process ensures that these surfaces are sufficiently alike that 
it is sensible to perform averages over results obtained from each individual surface. 
This similarity is summarized in the notion that the sample surfaces all belong to a sin- 
gle statistical ensemble. If measurements of the quantity f are made at the corre- 
sponding points (x,y) on each surface and the results are averaged to obtain a result 
f (x,y), there is reason to identify f (x,y) with the quantity 

(f(h,x,y)> = f f(h,x,y) p(h;x,y) dh (3) 

J h 

which is the expectation value of f(h,x,y) at the point (x,y). The probability density 
p(h;x,y) is that defined in equation (1) in terms of an idealized ensemble of sample sur- 
faces. Note that (f(h,x,y)) does not necessarily agree with any particular measurement 
at the point (x,y) on any particular experimental surface. A logical second set of experi- 
mental measurements can be performed by finding the set of values f(h,x,y) at a large 
number of points (x,y) on a single experimental surface, say the kth one, and subsequently 
averaging these results to obtain the quantity F^. Certain further stipulations about the 
character of the random process allow the identification of f ^ with a quantity calculated 
from the probability density. These stipulations are summarized in the notion that the 
random process is ergodic (ref. 2). An example of an ergodic process is the stationary 
Gaussian random process (with an integrable autocorrelation function) for which all the 
joint probability densities of equation (2) are independent of the location of the origin of 
the coordinates x and y. The second-order joint probability density for this process 
has the exponential form 1 indicated subsequently. When the random process is ergodic, 
one has 


f k = <f) = \ f(h,x,y) p(h;x,y) dh (4) 

J h 

Although the physical quantity f(h,x,y) may depend on position, the ergodic stipulation is 
such that the averages over the random process obliterate this dependence. Thus, the 
expectation value (f) is independent of position, and f ^ or <f> is a global property of 
the particular sample surface and does not necessarily correspond to any value f(h,x,y) 
measured at a particular point on the surface. A measure of the degree to which a single 
measurement at a single point may differ from f^ is the variance <Xf of the quantity 
f(h,x,y) which is given by 

Virtually all the work done to date pertaining to rough-surface scattering theory 
has assumed that the random processes involved are of this kind. 
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(5) 


a f 2 = <f 2 > . < f >2 

These statistical ideas enjoy successful application in many areas of theoretical 
physics, but their application to the problem of radar scattering from planetary surfaces 
is not without conceptual difficulties. The goal of many of the planetary scattering appli- 
cations (e.g., refs. 3 and 4) has been to study as small an area of the planetary surface 
as possible within the resolutional capabilities of the radar system, with the quantity <f> 
being related to the electromagnetic power backscattered by the area. In order for such 
an identification to be reliable, it is necessary first of all that the surface topography be 
describable as a sample function of a stationary ergodic random process. Unfortunately, 
terrestrial experience shows that natural terrain, in contrast to the ever-changing sur- 
face of the oceans, is fixed and over wide areas can assume large-scale features that 
obviously are not random in character. The fact that a measured quantity, such as back- 
scattered power, may fluctuate from point to point on such a surface is not evidence of 
randomness. An example which amply illustrates this is the rapid fluctuation with aspect 
angle of the radar cross section of an aircraft, a target whose every contour represents 
careful design rather than randomness. Thus, the basic assumption that a planetary sur- 
face is describable by a stationary ergodic random process is open to serious question, 
and certain contrary evidence is in fact in existence. For example, Hughes (ref. 5) has 
emphasized that a Gaussian height distribution is not appropriate for the Moon, and, as is 
discussed in the section "Experimental Results for Radar Scattering," it has been shown 
(ref. 6) that certain lunar radar studies have led to absurd results for the surface corre- 
lation length, a crucial parameter in the statistical description of any surface. 

An additional weakness in the statistical analyses surrounding planetary radar work 
is that apparently no one has calculated the variance of the scattered power (eq. (5)) in 
order to obtain an idea of the theoretical error associated with the reported results. 

The next section outlines the electromagnetic theory into which these statistical 
ideas are inserted to yield the overall planetary scattering theory. 

ELECTROMAGNETIC SCATTERING APPROXIMATIONS 

If S represents a closed surface containing all sources of the electromagnetic 
field, then it can be shown (ref. 7) that the electric field E p at a point P exterior to 
S is given by 


E 


P 


— © -j a>jLt(n x H)\j/ + (n x E) x Vi p + (n • E)Vm dS 
Jg L J 


( 6 ) 
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where E and H are the electric and magnetic field vectors on the surface and ft is 
the outward unit vector normal to S. Green's function = exp (-jk • r )/r, where r 
is the distance from a given point on the surface to the field point P. All fields are 
assumed to have an exponential time dependence exp (ju>r). Equation (6) can be trans- 
formed (ref. 7) into an alternative form 

Ep = --jL £ U |I - E dS (7) 

* 4tt Jg \ 3n 3n j 

Equations (6) and (7) are the mathematical representation of the Huygens -Fresnel princi- 
ple, which states that the field at a given point is produced by the sum of the fields of 
individual wavelets propagating from each point on the preceding wave front. Equation (7) 
can be written separately for each component of E, and the resulting scalar equation is 
the basis for much of the work in the field of physical optics. 

In order to apply equation (6) or (7) to rough -surface scattering, it is necessary to 
identify the closed surface S and to relate the fields E and H on S to their 
sources. Usually, S is identified with the actual scattering surface, and the source of 
the fields E and H is taken to be an isolated point, exterior to S, but placed so that 
no energy from the source can reach P without being scattered by S. Even if the 
shape of the scattering surface is known precisely, the determination of E and H on 
the surface in the presence of electromagnetic radiation from a distant point source is a 
formidable problem. The problem is greatly simplified if the scattering surface is a 
perfect conductor whose surface undulations are gently varying on a lateral scale of very 
many wavelengths. In such a circumstance the region in the neighborhood of any given 
point has the characteristics of a very large, nearly planar surface. Provided that there 
are no paths by which energy reflected from a given point on the surface can reach 
another point on the surface, the electromagnetic fields at the given point are approxi- 
mately those which result from reflection of a plane electromagnetic wave by a perfectly 
conducting plane tangent to the surface at the given point. When this approximation, 
called the tangent -plane approximation, is incorporated into equation (6) or (7), it forms 
the basis for most so-called physical -optics methods of surface scattering study. 2 For a 
surface satisfying the assumed conditions, the E and H fields in this approximation 
are easily evaluated from the incident field values and the boundary conditions for a per- 
fect conductor. Thus 


2 Rea, Hetherington, and Mifflin (ref. 8) have even asserted that the tangent -plane 
approximation is universally valid for perfectly conducting surfaces whatever their shape. 
However, a surface shape allowing a previously reflected ray as well as the incident plane 
wave to reach a given point will alter the total incident field and thus invalidate the 
approximation. 
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( 8 ) 


E • n = 2Ei • n 
E x n = 0 
H • n = 0 
H x n = 2Hi x n 

where Ej and H[ are the components of the incident plane wave. 

The physical-optics method with incorporation of the tangent -plane approximation 
is thus capable of describing the scattered electric and magnetic fields of a relatively 
smooth, perfectly conducting surface in a way that retains the coherence effects due to 
phase interference among field amplitudes scattered from different points on the surface. 
Under appropriate conditions^ which are rigorously discussed in references 7 and 9, the 
physical -optics formalism reduces to that of geometrical optics. Geometrical -optics 
propagation is characterized by the fact that the wave amplitude and phase at a given 
point on the wave front are related uniquely to those points and only those points on pre- 
ceding wave fronts which lie on the ray trajectory traced by the Poynting vector passing 
through the given point. The method of stationary phase, when applied to a wave front 
satisfying the conditions of geometrical -optics propagation, then leads to the results of 
ray optics (ref. 7) for which all diffraction phenomena are excluded. 

A theory to describe scattering from rough, natural surfaces could, on the basis of 
the present discussion, follow one of four procedures. First, equation (6) or (7) could be 
used, so that the results conform to exact methods of electromagnetic theory. For rough 
surfaces, this approach is entirely intractable. 

In the second approach the tangent -plane approximation for evaluation of the total 
surface fields could be applied to equation (6) or (7) so that the results contain the effects 
of point-to-point coherence. Unfortunately, most natural surfaces are dielectric rather 
than perfectly conductive. Even for plane waves incident upon locally plane regions, the 
reflection coefficients are complicated functions of the angle of incidence and are differ- 
ent for the two components of polarization. Furthermore, since a part of the incident 
wave energy penetrates a dielectric material, one must stipulate that none of this pene- 
trating energy can be scattered from underlying strata to emerge subsequently from the 
surface. 

The third approach would involve incorporation of the tangent -plane approximation, 
required assumptions for valid application to dielectric surfaces (as in the second 

^These conditions are satisfied in the far -field zone of a source or scattering object 
for which case the field variations both along and transverse to the wave front are small 
on a scale of many wavelengths. 
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approach), and in addition, application of a method of stationary phase, so that the results 
take on the character of a geometrical -optics approximation. This is the approach taken, 
for example, in references 10 and 11 and applied in references 3 and 4. Kodis in refer- 
ence 12 applies this approach in a consistent way for the special case of a perfectly con- 
ducting surface and effects the transition from equation (6) through the classical method 
of stationary phase to the formalism of pure ray-optics scattering. 

The fourth approach would invoke ray optics from the outset and would treat the 
scattering surface as a collection of mirrorlike facets from which the incident rays are 
specularly reflected according to the usual laws of elementary geometrical optics. This 
approach is followed in the work of references 8, 13, and 14. 

Each approach requires that the rough surface be treated in a suitable statistical 
manner. Details of the statistics vary among the approaches and give rise to certain 
difficulties in reconciling the results of the several approaches. First, that approach 
containing the least number of theoretically objectionable features, the method of Kodis 
(ref. 12) and certain extensions of his work by Barrick (ref. 15), are examined. 


STATIONARY -PHASE APPROACH OF KODIS 


Kodis (ref. 12) begins his work on perfectly conducting, rough surfaces by applying 
the tangent -plane approximation to an expression similar to equation (6) and by evaluating 
the surface fields by means of the equivalent of equation (8). He considers the incident 
plane wave to have unit amplitude when it reaches the surface, time dependence jwr, 
propagation vector k\, and electric field polarization direction e. He finds for the 
field E sc scattered in the direction of k s the expression 


® SC (ki,k s ) = exp (jkr) f F(n,ki,k s ) exp j(k\ - k s ) • r ’ 


dS' 


(9) 


where the integration is over the scattering surface, r ' is the vector position of a point 
on the surface, r = |r | represents the distance from the origin to the field point, and 

F(n,ki,k s ) = (T - kikg) • (iq x e x ft) 

with I being the unit dyad. For backscatter, ET S = -iTj, and F reduces to (ft • kjje. 
For sufficiently small wavelengths, the phase of the exponential term of the integrand in 
equation (9) varies very rapidly on a distance scale for which the remainder of the inte- 
grand varies only slightly. The principal part of the integral of equation (9) can thus be 
approximated by the classical method of stationary phase (refs. 7 and 16). Because of the 
rapid oscillation of the integrand, significant contributions to the integral can occur only 
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at the end points of the interval or at those interior points at which the derivative of the 
exponential argument vanishes. The exponential argument can thus be represented by a 
Taylor's series in the neighborhood of these stationary points, and the slowly varying 
remainder of the integrand can be evaluated at the stationary points themselves and thus 
removed from the integrand. The integral then becomes a summation of contributions 
from each of the stationary points. 

Kodis treats in detail the case of a one -dimensionally rough or corrugated surface 
with height in the z -direction given by h(x) (£ in Kodis' notation) and width L in the 
y-direction extending from y = -L/2 to L/2. He also gives the results of the straight- 
forward extension to a two -dimensionally rough surface. In both cases the mean surface 
is assumed to be planar. He finds for the one -dimensional case 


N 

E sc (k i ,k s ) ^ 3!L - Fm § exp (jkr) J exp (jkq . r ’) dy' dx' 
m = 1 m z m 


N 

y i k F rn 

i 2™ ft m • g z 
m=i 


exp 



exp jk 


? in 



- x m ) 2 • r’)| x , =x , 

dx’ 2 x m 


dy’ dx’ 


( 10 ) 


where q = k\ - it s , dS' has been replaced by dy' dx'/(n m • e z ), and the points x^ 
are the points at which q • F is stationary with respect to x and y. Thus, at these 
points 


dh(x') 


dx’ 


l m 




(ID 


Equation (11) implies that q is collinear with n m , the unit normal at the stationary - 
phase points x^. This is the usual condition for specular reflection. The integrals in 
equation (10) are performed over small neighborhoods S m about each of the points r^ 
and can be written as 
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I m = exp (jkq • rin)^ exp jk 

s m 


q z h"(xm) 


dy’ dx' 


= I-> exp (jkq • r l exp jk 
JAx m 


(x - x{n ) 2 


h "(xm) 


dx T > 


= L exp (jkq • r Jn) \ exp jk 
^Au 


Y q z h " (x m) 


du 


where u = x' - x^, Au - Ax m is a small interval around x^, and 


(12) 


h"( x m) = '^~9* 1 (x') 
dx’ 2 


x’^xJn 


In the last integral in equation (12), the increasingly rapid oscillation of the integral with 
increasing u renders the integral insensitive to the exact interval Au, so that the 
stationary-phase approximation allows Au to be replaced by the interval (- 00 , 00 ). Then 
use of the formula from reference 16 


p°° 

j exp (jm2u2) du = exp (j ^ 


reduces equation (12) to 


I m = L exp Jj exp (jkq • r' m ) 


2tt 


kq h”(x^) 


1/2 


(13) 


Substituting equation (13) into equation (10) gives for the absolute square of the scattered 
electric field the value 



2 


= E sc - 





where use has been made of the fact that at the specular points n m = -2.. The total 
radar cross section 2 is then related to the scattered field by the expression 


14 



S(kj,k s ) = 47rr 2 | E sc | 2 


and the cross section a per unit_width L of the corrugated surface is then 


o(?i,E- s ) = L^|F 

a 0 
H z 



h"<x|n) h”(xj,) 


- 1/2 


exp 



( 14 ) 


where the ensemble average has been inserted and is taken over the physically possible 
realizations of random surfaces. The number of specular points N, the positions r ^ 
of the specular points, and the quantities h'^x^) evaluated at those points differ for 
each realization. 

It is straightforward to show that for the surface under discussion, 



(15) 


where p m is the radius of curvature at the point x^. Also, use of equation (11) gives 



(16) 


Substitution of equation (16) into equation (15) and of that result into equation (14) gives 


N 

£ {PmPp ) 1/2 ex P 

P=1 
p*m 

(17) 

The first ensemble average on the right side of equation (17) can be written 





(18) 
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where p m is the average radius of curvature at the N specular points of a single 
realization of the ensemble of rough surfaces. Kodis makes the assumption that although 
the number N of specular points on a given surface realization may vary from one 
realization to another in the ensemble, the average radius of curvature at the specular 
points p m is the same for all realizations of the ensemble. Since the angle of inci- 
dence of the radiation, and hence the locations of the specular points, is arbitrary, this 
assumption is equivalent to assuming that the random process generating the radius of 
curvature at any point on the surface is an ergodic process. If this assumption is true, 
then 


(Np m > = <N>P m 


(19) 


An equation analogous to equation (19) can be written for the second ensemble average of 
equation (17) as 




<N(N-l)>(p m p p) l/2 exp 



(m * p) (20) 


Kodis uses equations (18), (19), and (20) in equation (17). The physically plausible 
assumption that the location of a specular point and the radius of curvature at that point 
are independent random variables over each surface realization allows Kodis to write 


(PmPp ) 1/2 ex P i (*m - *p) = {PmPp ) 1/2 ex P 3 (<t> m ~ *p) ( m * P) (21) 


where (f>^ = kq • r { with t = m and p. For a large number of specular points N and 
with <p m - <p p distributed uniform ly over all angl es (as is the case for any given surface 
if the wavelength is short enough), exp j(4> m - 0 p ) on the right side of equation (21) 
tends to zero. Thus, equation (17) reduces to 

cr(ki,k s ) = L ^|F| 2 (N>p m (22) 

For the backscattering case k s = -kj, |F| = 1, and q = 2. Kodis extends his calcula- 
tions from one -dimensionally rough to two-dimensionally rough surfaces and obtains for 
the total backscattering radar cross section 
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S(fi,* s > = N> 


(23) 


where (N) is the ensemble average of the total number of specular points and p a p^ is 
the average value over a surface realization of the product of the two principal radii of 
curvature at each specular point. 

Equation (23) is analogous to a result given In reference 7 for scattering from a 
smoothly curving object in the stationary -phase or geometrical-optics approximation. It 
is also analogous to the geometrical -optics cross section of a smoothly curving body, 
such as an ellipsoid or sphere, given in reference 17. The theory leading to equation (23) 
predicts no depolarization of the reflected wave, so that if depolarization of the signal is 
actually observed, there is reason to believe that the underlying assumptions for validity 
of the theory are invalid. Depolarization in actual planetary radar reflections is dis- 
cussed further in a subsequent section. 

BARRICICS EXTENSION OF KODIS' RESULTS 


In order to complete the calculation of the radar cross section (eq. (23)) it is nec- 
essary to evaluate the quantities (N) and Some work of Barrick (ref. 15) that 

attacks this problem directly is discussed here. Barrick notes that the problem of 
counting the number of points on a corrugated surface at which the slope attains a speci- 
fied value (i.e., the number of specular points on the surface) is formally similar to the 
problem of counting the average number of zero crossings per unit time of a random 
function of time with specified statistical properties. This problem is solved by 
S. O. Rice (ref. 1), and its solution is adapted by Barrick to rough surfaces and extended 
to a two-dimensionally rough surface. For the two-dimensional surface, Barrick finds 
for the expected number n^ = <N> /Total surface area of specular points per unit area 


p°o /-»°o poo ! 2 I 

n A = P( h x spjhy sp) j ^ j m j w I ^xx^yy - h x y|p^h xx ,h x y,hyy|h x sp,hy spj ^xx ^xy ^yy 

(24) 

where 


h, - 9h 
% sp - 9X 


x=x 


m 


h = 9h 

n ysp 9y 


y=yln 


and where p(h x sp>hy sp) is the probability density for the surface slopes evaluated at 
the specular values and p^hxxjhxyjhyyjhxgpjhygpj is the conditional probability density 
for the three second derivatives of the surface height given that the first derivatives have 


17 



their specular values. The double subscripted quantities hxx, h xy , and h yy corre- 
spond to the respective second partial derivatives of the surface height h(x,y). The 
functional notation for h has been dropped for convenience. 


Barrick then develops an expression for p a p^, or ^|rjr 2 |^ in his notation. (This 
notation implies that he is replacing the average over a single surface with an ensemble 
average, a permissible substitution under a suitable ergodic hypothesis.) In his develop- 
ment, he adopts the rather unsatisfactory tactic of replacing with — 1 - where 


|J| = 


hxxhyy " hxy 


<UI> 


Barrick asserts that no quantitative measure of the error involved 
in this replacement is available, but that for large values of | j|, the error should not be 
excessive. It is shown here by an alternative method that Barrick' s final result for 
(| r l r 2|) exac *> so that his unsatisfactory assumptions regarding are not 

essential. 


From equation (24), the ensemble -averaged number dn^ of specular points per 
unit area such that at the specular points the second derivatives of the surface have 
values within dh^ of h^, within dh xy of h xy , and within dh yy of h yy , is 


dnA - p(hxsp>h y sp) hjochyy - hxy P^h xx ,h xy , h yy |h x S p,h y S pj dhxx 


dh X y dhyy 


(25) 


Assume that dn^ also represents the number of specular points with the above pre- 
scribed values of h^, h xy , and h yy for a single surface realization from the 
ensemble. This will be a plausible ergodic assumption if the number dn^ of such 
specular points is large. At these specular points the quantity j r ir 2 | can be written 


I r l r 2| = 


(l + hx sp + hy sp) 


hxxhyy 


xyl 


sec** y 

hxx^yy - h xy 


(26) 


where y is the angle between the mean surface normal and the local surface normal. 
To find the average of |r jr 2 | over a single surface, it is necessary merely to add up 
the quantities represented by equation (26), each quantity being multiplied by the number 
of times it occurs (from eq. (25)), and then to divide by the total number (from eq. (24)). 
Thus 


r l r 2 


J | r l r 2| dn A 
"A 


sec^ y yyy p(h X x>hxy>hy y | hx 

spjhysp) d ^*xx d h xy dhyy 

yyyjhxxhyy - h X y| P^hxxjhxyjhyy 

hx sp>hy spj d hxx dh xy dh yy 


(27) 
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Since the triple integral in the numerator of the right member of equation (27) is 
unity, equation (27) is identical with the final result of Barrick. Dividing equation (23) by 
the total surface area and substituting equations (24) and (27) give for the radar cross 
section per unit area 

a = n sec 4 y p(hx sp ,h y sp ) = w sec 4 y p(y) (28) 

where p(y) = p(h x S p,hy sp j . At each specular point contributing to the result in equa- 
tion (28), the direction of incidence is parallel to the local surface normal. It is plausi- 
ble therefore that equation (28) can be extended from a perfectly conducting surface to a 
dielectric surface by insertion of the reflection coefficient for normal incidence on a 
dielectric plane. Note that such a procedure is justifiable only after the reduction of the 
problem to that of a summation of scattering contributions from the set of specular points. 
An attempt to insert into equation (9) a dielectric -surface reflection coefficient, which is 
a complicated function of the local angle of incidence at each point on the surface, would 
have resulted in additional complicated factors in the integrals over S m of equation (10). 
These factors would have altered the subsequent development substantially and would have 
led to problems of interpretation similar to those discussed in the section "Correlation- 
Function Method." 

Barrick quite plausibly inserts the dielectric reflection coefficient into equation (28) 

to get 


7 tv 


= t t sec 


y p(y) 



(29) 


as the backscattering cross section per unit area for a dielectric body, where £ and rj 
refer to the polarization states (horizontal or vertical) of the incident and scattered waves, 
respectively, and the absolute values of the Fresnel reflection coefficients for normal 
incidence of a plane wave on a plane surface are denoted by |R^|. These coefficients 
are such that 



where h and v refer to horizontal or vertical polarization. Equations (29) and (30) 
show that for a dielectric as well as a perfect conductor, there is no depolarization of the 
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radar signal upon reflection in this approximation. If depolarization is actually observed, 
these formulas must be used with extreme caution. 

Equation (28) - which agrees with the results of Muhleman (ref. 13), who based his 
work from the outset on a pure geometrical -optics (ray-optics) method - has thus been 
traced back to first principles, so that its domain of validity is reasonably clear. 

Although equation (28) or (29) indicates that the angular scattering dependence is the 
same as the distribution of surface slopes, this may not be true for arbitrary slope dis- 
tributions. In particular, the slope distribution must fall off with y sufficiently rapidly 
so that neither surface self -shadowing nor multiple scattering is important. Surface 
self -shadowing would tend to reduce the number of contributing specular points per unit 
area and would be important for large angles of incidence or large surface slopes. Sim- 
ple geometrical considerations show that multiple scattering could be important if there 
are a significant number of surface slope angles a such that a ^ | For slopes 
with such angles, incident rays can reach a given specular point by intermediate reflec- 
tion from nonspecular points located on those slopes. 

An alternate route to equation (29) has occupied a principal place in the literature 
of rough-surface scattering. That route is through the so-called correlation -function 
method, whose intermediate results have provided many of the formulas used in planetary 
radar data reduction methods. 

CORRELATION -FUNCTION METHOD 

The theoretical work which has had the most impact on radar studies of planetary 
surface properties has probably been that of Hagfors (refs. 10, 11, and 14). This work, 
which falls into the third category discussed previously in the section "Electromagnetic 
Scattering Approximation," unfortunately contains many contradictory features, and appli- 
cations of parts of it to some lunar and planetary radar data analyses (refs. 3, 4, and 18) 
have led to unreasonable results. Some of these results are discussed in the following 
section. The purely theoretical aspects of Hagfors' work are discussed in this section. 

In reference 10, Hagfors attacks the problem of plane -wave scattering from a one- 
dimensionally rough (corrugated) surface with random height deviations from a mean 
plane. He extends the final result to the two-dimensionally rough, mean-planar surface. 
His work for this special surface is presented in detail, since it will help to clarify his 
work (ref. 11) on a rough spherical surface which is the basis for certain data analysis 
techniques. In Hagfors' development, the height of the rough surface in the z -direction 
h is independent of the y -direction which extends from to +°°. The slope of the 
surface in the x -direction is denoted by t, so that the slope angle a is given by 

a = tan - * t. The angle of incidence of the impinging wave with respect to the mean sur- 
face is denoted by cp. The local angle of incidence is then given by \p - (p - a. As 
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before, k is the wave number 2 n/\. The scalar value of an electric field component of 
the incident wave at point x on the surface is then of the form 

E = E 0 exp [-jk(x sin cp - h cos <p)] (31) 

By means of the tangent -plane approximation with the necessary stipulations for applica- 
tion to dielectric bodies, the scattered electric field amplitude on the (dielectric) surface 
is given by E Q(i //), where Q(i//) is the plane-wave reflection coefficient. After 
straightforward manipulation of the scalar form of equation (7), Hagfors writes for the 
backscattered field amplitude E R at the radar antenna 


E R = exp (-jkR) y Q(^) cos ^ exp [-2jk(x sin q> - h cos <?)] ds (32) 

where R is the distance from the origin to the radar antenna and ds is the line ele- 
ment along the rough surface. Introducing dx = cos a ds and eliminating i]/ in favor 
of t and (p , Hagfors writes 

E r = K f Q(t)(cos (p + t sin <p) exp £-2jk(x sin (p - h cos <p)j dx (33) 

where K represents the constant factor outside the integral in equation (32). The power 
received at the radar antenna is then proportional to and the expectation value of 

this quantity over a suitable ensemble of surfaces describing the random height h is 
proportional to the ensemble -averaged received power. Thus 


< E R E r) = (|%| 2 ) = |K|2 II <*1 dx 2<^( t l) Q*(t 2 ) (cos <P + ti sin <p) (cos (p + t2 sin <p) 


exp 2jk(hj - h 2 j cos <p\ \ exp -2jk(xj - X 2 ) sin <p 


where the subscripts 1 and 2 denote dummy variables of integration. Introduction of the 
new variables Ah = hj - h 2 and Ax = xj - X 2 gives 
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( 35 ) 


(|e| 2 ) = |k| 2 j* dxi j* d(Ax)^Q(ti) Q*(t 2 ) (cos cp + 1 1 sin <p)(cos <p + t 2 sin <p ) 
exp (2jk Ah cos cpf^ exp (-2jk Ax sin <p)| 


In order to evaluate the expectation value within the integral in equation (35), a joint 
probability density function corresponding to equation (2) must be introduced for the 
height difference Ah and the slopes tj and t 2 . Hagfors assumes this joint probabil- 
ity density to be that of a Gaussian process and writes explicitly 


p(Ah,tit 2 ) = (27r) _3 / 2 (det A)"V 2 





+ 2M 12 Ah(t x - t 2 ) + 2M 23 t 1 t 2 } 




(36) 


where the Mjj are the cofactors of the symmetric matrix 



2[l - p(Ax)] 

9p(Ax) 
9 (Ax) 

9p(Ax) 
9 (Ax) 


9p(Ax) 

9 2 p(Ax) 


9 2 p(Ax) 

2 

l rms 

8 (Ax) 

9 (Ax) 2 

Ax=0 

9 (Ax) 

I 2 


9p(Ax) 

9 2 p(Ax) 

9 2 p(Ax) 



9 (Ax) 

9 (Ax) 2 

9 (Ax) 2 

Ax=0 


(37) 


and where det A is the determinant of A. The quantity p( Ax) is the autocorrelation 
function of the random surface height and hp ms is the root-mean-square height. Ref- 
erence 1 or 2 cover the theory surrounding equations (36) and (37). To use this formal- 
ism, it would be necessary to identify the expectation value in the integral of equation (35) 
with the quantity <f> of equation (3) and to substitute equation (36) into equation (3), where 
the single integral had been replaced by a triple integral over Ah, t j_, and t 2 . 

At this point a hard fact comes to the surface, namely, that Q(t) is a very com- 
plicated function of t which has a different functional form for each of the two possible 
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local polarizations of the incident wave. Since the local polarization changes with posi- 
tion according to the detailed variation of the rough surface, it is completely unclear how 
to proceed from this point. Hagfors' approach is to seek a method for effectively remov- 
ing the t dependence of Q. His argument is as follows: Suppose that certain slope 
values ti and F 2 can be found such that the quantity F(t^,t 2 ) defined by 


F ( t l» t 2) = Q(*l) Q*(*2) ( cos ^ + *1 sin <P) ( cos <P + *2 sin <?) (38) 

can be expanded in a Taylor's series about the points Fj and t 2 with the first deriva- 
tive terms vanishing identically. This condition is met if 


y dxj y d(Ax) exp (-2jk Ax sin (p) d(Ah) dtj dtj 


(*i 


ti~ti 


- ti) exp (2jk Ah cos cp) p(Ah,ti,tj) 


= 0 


(39) 


where i,j = 1,2 with i * j and p is defined by equation (36). 

Hagfors asserts that equation (39) yields 

ti = t 2 = tan <p (40) 


If equation (40) holds, and if the second and higher derivatives are arbitrarily dropped 
from the Taylor's series for F, then F reduces to 

F ( t l> t 2) = F (tl>*2) = Q( tan •?) Q*(tan cp) cos -2 (p (41) 

The value Q(tan cp) then corresponds exactly to the Fresnel reflection coefficient for 
plane-wave normal incidence on a plane surface, which is independent of polarization. 
Thus, Q(tan <p) Q*(tan <p) is given by equation (30) (with j u. set equal to unity). Use of 
equations (30), (40), (36), and (37) in equation (35) then gives 


(| e r| 2 ) = |K| 2 ^ ■ ^ J (fai j* d(Ax) exp (-2jk Ax sin <p) j* dtj J dt 2 j* d(Ah)|exp (2jk Ah cos <p) p(Ah,ti,t 2 )j (42) 
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Equation (42) can now be integrated over tj, t 2 , and Ah to give 


/|Er| 2N > = \k\2 (— — f — - — f dx x f d(Ax) exp (-2jk Ax sin <p) 

\|R|/ 1 \Ve + l) cos 2 <p J -L/2 1J - 

exp jj4k 2 hp ms cos 2 (p \_ 1 - p(Ax)]J 


( 43 ) 


where L is the extent of the surface in the x -direction. Equation (43) can be extended 
to a two-dimensionally rough surface for which p(Ax) can be identified with p(Ar), 
where Ar is the separation distance between two points on the surface. The one - 
dimensional Fourier transform indicated by the integration over Ax in equation (43) 
can be replaced by a two-dimensional Fourier transform such that 




exp 


-4k 2 h 2 ms cos 2 <p[l - p(Ar)| J 0 (2k Ar sin cp) 


(44) 


What Hagfors' method actually represents is a form of stationary -phase approxi- 
mation similar to Kodis' method but different in detail from this classical method. The 
previous discussion shows that the particular choice t[ = tan (p has nothing to do with 
the ensemble average indicated in equations (34) or (35). If the ensemble -average brack- 
ets of equation (34) are moved outside the integrals, the two exponentials of the integrand 


may be combined to give -r^ = = tj = t2 = tan <p as a point of stationary phase of the 


dh x _ dh 
dxj dx2 

integrand. The method of stationary phase permits the remainder of the integrand to be 


evaluated at this point and removed from the integrals. Up to this point, the procedure 
is the same as that of Kodis. Instead of evaluating the remaining integrals only in the 
neighborhood of the points of stationary phase, however, Hagfors* method retains the 
entire integration range. When the ensemble -average brackets which were placed out- 
side the integrals of equation (34) are moved back within the integrals over xj and Ax, 
the result is equation (42). 


This view of equation (43) shows it to be a combination of the method of Kodis and of 
an attempt to retain the effects of point-to-point coherence of the scattering process. The 
coherence retention is not complete, however, since the reflection coefficient allows only 
for locally normal incidence. Therefore, coherence can be retained, at best, only among 
the isolated specular points, but the integral on Ax of equation (43) contains the phase 
variations of a continuum of points over the entire surface. The physical meaning of 
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equation (43) is thus somewhat obscure and it is unclear whether it is capable of describ- 
ing to any reliable extent the effects of surface correlation or coherence. 

ANGULAR SCATTERING LAWS 

In order to proceed from equation (44) to a practical result for real surfaces, a 
form for the autocorrelation function p( Ar) must be assumed. For reasons of mathe- 
matical simplicity, attention in the literature has been devoted almost exclusively to two 
forms: the exponential correlation function 

p(Ar) = exp (45) 


and the Gaussian correlation function 


p(Ar) = exp 



(46) 


where l is the correlation length of the statistical process representing the rough sur- 
face. Actually, equations (45) and (46) are almost never used as they are written. For 
large values of kh rms cos <p, the so-called deep -phase modulation case, the exponential 
factor containing p(Ar) in equation (44) drops rapidly to zero with the slightest depar- 
ture of p(Ar) from unity. It is generally assumed that equations (45) and (46) can be 
approximated by the first two terms of a Taylor's series expansion about a point at or 
near the origin. A controversy has appeared in the literature over whether the point of 
expansion should be the origin or another point near it and whether an expansion of equa- 
tion (45), which has no Taylor's expansion about the origin, can represent a sensible cor- 
relation function at all. This controversy is discussed in detail by Barrick in reference 6 
and to some extent by Fung and Chan in reference 19. 

If equation (45) is approximated by p(Ar) = 1 - Ar/j? and if this form is substituted 
into equation (44), the latter equation can be integrated exactly to give for the power 
reflected per unit surface area 

2 I " 3 / 2 

cos^ <p + - — •) sin 2 (p (47) 

(2kh rms ) 2 \ h rms/ 


<P R (<P)) 


v/£ - i 

x/e + 1 
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If equation (46) is approximated by p(Ar) = 1 - (Ar) 2 /4 2 , the analogous result is 

<P R <*)> a M. - \ — exp ^ — tan 2 (48) 

\ R ; U + V COS 4 CP \ 4h 2 ms ) 

Hagfors (ref. 10) obtains these results for a mean -planar surface and in reference 11 he 
ostensibly obtains the same results for the case for which the mean surface is sperical 
as is needed in planetary studies. The development of reference 11 is somewhat unsat- 
isfactory, however, since Hagfors is forced to drop the distinction between the sub- 
scripted dummy variables in the analog of equation (34). Furthermore, he makes some 
approximations in his integrand functions that are valid only if <p = 0. The effect of 
these arbitrary steps is to suppress the spherical nature of the surface. Notwithstanding 
these incongruities, he asserts that equations (47) and (48) represent the expected varia- 
tion with incidence angle of the power backscattered from a spherical planetary surface, 
and these expressions are used in actual data reduction of radar scattering from Mars in 
references 3 and 4. 

Note that the use of either equation (45) or (46) in equation (44) produces an integral 
over Ar that does not exist because of the infinite upper limit of that integral. It is 
only the two-term approximations which lead to a finite result. Fung and Chan (ref. 19) 
obtained numerical results for the integral in equation (44) by using an unspecified upper 
cutoff D for that integral. Although they did not comment on the dependence of their 
result on the value of D, it can be readily shown that their numerical results do, in fact, 
depend on D. Their plotted curves result from a specifically chosen value of D that 
bears no apparent relationship to any other physical parameter of the problem, and their 
results would have diverged as \/D for larger values of D. 

Extensive difficulties with the exponential correlation function of equation (45) have 
been discussed by Barrick (ref. 6). He points out that this surface correlation function 
violates the fundamental assumption of the theory that the surface is of a gently varying 
nature. Equation (45) corresponds to a surface with many abrupt, vertical discontinuities. 
Furthermore, Barrick shows that in the limit of small radar wavelength, which is the 
only limit in which the developed theory can claim validity, the total power scattered in 
all directions by a perfectly conducting, rough surface described by equation (45) is zero, 
even if the surface slopes are small. Equation (46), on the other hand, describes an 
admissibly smooth surface and leads to a plausible value for the scattered power, if the 
slopes are small. 

The angular scattering law of equation (48), being based on equation (46), does not 
lead to the absurdities pointed out by Barrick, although it is still subject to the inconsis- 
tencies discussed in the previous section. As e — °°, the dielectric -surface reflection 
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coefficients for all angles of incidence reduce to those appropriate foy a perfectly con- 
ducting surface. Similarly, if the limit e - «° is taken in equation (48), the resulting 
angular scattering law agrees with the angular scattering law deduced by Beckmann 
(ref. 20) for a perfectly conducting, rough surface having the same statistical properties. 
Although no satisfactory derivation of equation (48) for a dielectric surface has ever been 
given, the agreement of that equation as e — 00 with a correctly derived result for the 
perfectly conducting surface makes it conceivable that there does exist some undiscovered 
basis for its validity. From the standpoint of Martian radar studies, however, this con- 
jecture is merely academic, since the experimental results are not describable by 
equation (48). 

EXPERIMENTAL RESULTS FOR MARTIAN SURFACE SCATTERING 

The Haystack Microwave Facility of the Massachusetts Institute of Technology 
Lincoln Laboratory has been extensively used to study radar signals reflected from the 
surface of Mars (refs. 3 and 4). In these studies, determination of the dielectric constant 
and surface slopes of Mars has relied upon the previously discussed theoretical work of 
Hagfors. It is unfortunate that, within the context of this theory, only the theoretically 
unacceptable correlation function of equation (45) has been found to fit the experimental 
data. 

The data reduction procedure used in references 3 and 4 consisted of fitting the 
angular scattering law of equation (47) to the experimental angular scattering function 
and matching the theoretical and the experimental powers at (p = 0 and at another value 
of <p . The fit at the two points gives a value for the dielectric constant e and an exper- 
imental number for the coefficient of sin 2 (fiA This coefficient, denoted by C, is such 
that C -1 / 2 i s approximately equal to the tangent of the angle (p for which equation (47) 
falls to 1/e of its value at cp = 0. Although this is called the root-mean-square sur- 
face slope in reference 4, this designation is clearly a misnomer, since for a surface 
described by equation (45), which contains arbitrarily large slopes with many vertical 
discontinuities (ref. 6), the root -mean -square slope is infinite. Arguments for the inter- 
pretation of as an "effective" slope are given in references 3, 11, and 14. These 

arguments center around the assertion (ref. 14) that the effective slope seen by the radar 
depends upon "... structural detail in the correlation function near the origin which will 
not appreciably influence the value of the integral [equation (44^| determining the back- 
scattered power." Unfortunately, it is precisely the form of the correlation function near 
the origin which gives rise to the angular scattering law used (eq. (47)). A change in the 
behavior near the origin from a linear variation in Ar to a quadratic variation is just 

4 In addition to radar system parameters, the constants of proportionality of equa- 
tions (47) and (48) also contain the unknown quantities ft and h rms , so that e cannot 
be obtained by a single measurement at cp =* 0. 
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the change necessary to achieve equation (48), which does not fit the data, rather than 
equation (47). Ignoring this point, Hagfors argues that the finite radar wavelength 
imposes a "filtering action" on the wave number spectrum generating the rough surface, 
so that the higher spatial frequencies of this spectrum are truncated or attenuated. Thus 
the infinite actual slopes of the model are of no consequence and interpretation of C -1 /2 
as the "effective" root -mean -square surface slope is justified. This argument is also 
unsatisfactory, however, since such a filtered spatial frequency spectrum does not give 
rise in the first place to the angular scattering law containing the C parameter. 

Because of these basic theoretical difficulties, the reliability of the Martian 
dielectric -constant and slope values quoted in references 3 and 4 is unclear. It may be 
thought that an indicator of the reliability of the method for Mars might be application of 
the same method to lunar radar data. Such is not the case however, because of the fol- 
lowing considerations. Earth-based lunar data have given results for the dielectric con- 
stant in the range from 2.5 to 3, which agree well with actual values. These good values 
for the Moon are clouded, however, by other difficulties. Barrick (ref. 6) has pointed out 
that the lunar radar data imply the absurd result that the surface correlation length £ is 
as large as one -half of the distance from the Earth to the Moon. Notwithstanding this 
intrinsic contradiction, one might argue that a theory can give a right answer for the 
wrong reason, so that some confidence in the dielectric -constant values for Mars could 
still be retained by analogy with the lunar case. Such an argument could have merit only 
if the Martian dielectric constants and surface roughness were qualitatively similar to 
those of the Moon. The spread in radar -obtained values of e for Mars ranges from 1.5 
to more than 6, and the radar data indicate that Mars is ostensibly much smoother than 
the Moon. This marked qualitative difference in the radar return from the two bodies 
severely weakens the probability that a defective method, which happens to be successful 
for e for the Moon will be successful for e for Mars. 

Other radar studies of Mars, using the radar facility at the Goldstone Tracking 
Station, has been reported in reference 21. These studies produced radar cross-section 
data and an empirical angular scattering function, but since no a priori theoretical model 
was used, it was impossible to produce any local values of dielectric constant. It is 
notable that the empirical angular scattering function resulting from these studies does 
not agree with equation (47), the angular scattering function of references 3 and 4. 

DISCUSSION 

The difficulties of the correlation -function method tend to obscure the relationship 
of this method to other methods of analysis of rough -surface scattering. Hagfors has 
concluded (ref. 14) that for the Gaussian autocorrelation function, the correlation -function 
method is equivalent to that of equation (29). The status of the exponential correlation 
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function with regard to this equivalence remains unclear. To the extent that this equiva- 
lence holds for the Gaussian case, only one type of acceptable theory for radar scattering 
from rough, dielectric planetary surfaces, namely, the theory underlying equation (29), is 
left. This theory is constrained to small slopes and gently varying terrain on the scale 
of a radar wavelength. Therefore experimental data to which the theory could be reason- 
ably applied must necessarily fall off rather sharply with the incidence angle of the sig- 
nal (<p of eqs. (47) and (48) and y of eq. (29)). The Martian radar data fall off rather 
slowly with incidence angle, however, so that the theory of equation (29) is still 
inadequate. 

Refinements to the specular -point theory include the introduction of a shadowing 
function into equation (29) so that obscuration of specular points near the planetary limbs 
can be taken into account. Beckmann (ref. 22) has developed such shadowing functions 
and has applied them with some success to lunar scattering data, although Brockelman 
and Hagfors (ref. 23) contest the success of his method. In any event, the introduction of 
shadowing would tend to drive down the level of the return with increasing angle of inci- 
dence, an effect which is opposite to that desired for explaining the Martian data. Also, 
inclusion of shadowing does not alter the limitation that the basic tangent -plane - 
approximation theory, to which the shadowing is appended, cannot take into account mul- 
tiple scattering or small-scale structure. These factors would be more important for 
regions near the planetary limbs, since their contribution to the total scattered power 
would not decrease with increasing incidence angle to the extent that the specular -point 
contribution does. The presence of these factors on the planetary surface would be char- 
acterized by larger -than-expected backscatter levels and increasing depolarization of the 
backscattered signal in the data returned from the planetary limbs. Both these features 
are invariably present in the radar data, giving further evidence that any theory based 
upon the tangent -plane approximation is incapable of describing backs mattering for large 
incidence angles. 

Other theoretical models for planetary radar data analysis are composite surface 
models containing two or more correlation -length scales. Some of these methods (e.g., 
ref. 24) are straightforward extensions of the basic theory described in the present 
report, whereas others (refs. 25 and 26) attempt to combine a large-scale roughness the- 
ory based on these same methods with a small-scale result based on a perturbation solu- 
tion to an exactly soluble scattering problem. Recent work by Burrows (ref. 27) has 
indicated that the small -perturbation solution is, in fact, equivalent in many respects to 
the solution obtained from the tangent -plane approximation, if this approximation is arbi- 
trarily applied to a surface with structure much smaller than the radar wavelength. 

Since in this case the physical basis for the validity of the tangent -plane approximation is 
completely lacking, the meaning of Burrows' result is unclear, but it does cast some 
doubt on composite models using a small -perturbation approach. None of these more 
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elaborate approximations have yet influenced the interpretation of Martian radar data. 
More disheartening, however, is the fact that not one of .the theories available at present 
is likely to bring order out of the confusing and sometimes self -contradictory radar 
results. 

The selection of references which have been discussed in detail has been made in 
order to present work that has had considerable impact upon planetary radar studies. 

In general, the best theoretical work has been limited to perfectly conducting surfaces. 

For example, Hoffman (ref. 28) gave in 1955 an excellent treatment of scattering from a 
rough, perfectly conducting surface in the tangent -plane approximation. He obtained for- 
mal results for the scattered Poynting vector, including the cross -coupling of the polari- 
zation components (depolarization) produced by local variations in the surface normal. 
Probably because of the formal nature of his results, Hoffman's work has had no impact 
upon practical applications. An obvious extension of Hoffman's work in the direction of 
the theories discussed in this paper could be incorporation of the reflection coefficients 
for arbitrary incidence angles in order to cover the dielectric scattering case. Even if 
successful, such a course would no doubt produce complex formal results not readily 
amenable to data analysis. In addition, the fundamental difficulty (discussed previously 
in the section "Statistical Nature of Rough-Surface Scattering") remains that planetary 
surfaces may not be describable in terms of the theory of random processes, so that 
attempts to extract meaningful results from such theories may be hazardous at best. 

CONCLUDING REMARKS 

The theoretical basis for the most widely accepted and used models for describing 
planetary rough-surface scattering has been examined. These models are basically more 
capable of describing rough, perfectly conducting surfaces than the dielectric surfaces of 
real planets. In order to obtain models which can describe dielectrics, it has been nec- 
essary to ignore all points of the planetary surface for which the local angle of incidence 
of the radar wave is not zero. Such a procedure cannot be justified unless the scattering 
surface is gently undulating with no structure on the scale of a radar wavelength. The 
presence of depolarization of the returned radar signal near the planetary limbs may be 
evidence that multiple scattering or small-scale structure actually exists, so that the 
theory, which predicts no depolarization, may not even be applicable to the sub- Earth 
areas of the planet that show no significant depolarization. 

In addition to the difficulty in extending the theory to dielectric scatterers, there 
are numerous questions concerning the theoretical development of the statistical scatter- 
ing properties. The theoretical formulation underlying the most often quoted experimen- 
tal results for the planet Mars has been found to be internally inconsistent and to lead to 
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an absurd result for the total scattered power in the limit of small wavelength for which 
the theory has its only clear basis for validity. 

The presently available approaches to the theory are incapable of accounting for the 
slow falloff of scattered power with angle of incidence for the planet Mars, and introduc- 
ing shadowing effects into the theory produces a correction in the wrong direction. In 
view of the difficulties involved with the basic theory, it is not apparent that more elabo- 
rate methods, such as the use of composite roughness scales, would produce better 
results. 

It is concluded that a fundamental breakthrough in the theoretical description of 
scattering from rough terrain is required before dependable results can be expected for 
such detailed properties as dielectric constant and surface slopes when even the statisti- 
cal properties for the scattering body remain unknown. Such a breakthrough represents 
a formidable problem that is not likely to be solved in the near future. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va., 23665 
May 2, 1975 
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